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Navigating the Slopes of a Methane Frost-Capped Mountain on Pluto 
(A Simulation) 


1 Introduction 

This application note simulates the speed of a drone as it travels down a 
mountain slope covered with methane frost. Initially, the drone is at rest atop 
the slope. Once its brake is released, the drone begins to descend under the 
influence of atmospheric drag (assumed to be quadratic) and sliding friction. 
In section 2, a mathematical model is developed that describes the drone’s 
speed as a function of time. 


For comparison, analytical, modified Euler method, and analog computation 
solutions will be determined and displayed in table 1. 


2 Mathematical modeling 
The following mathematics will be a bit lengthy, so you may want to grab a 
snack, some iced tea (stirred, not shaken), a comfortable chair, and hang on 
tight! 
Starting with Newton’s Second Law of Motion, 

ma = Fnet = -drag - friction + slope component of force due to gravity 

ma = -kv? — umgcos(@) + mgsin(@) 
a= -k/m v? + [sin(@) - pcos(8)|g 

Letting a = dv/dt and g’ = [sin(®) - ucos(@) |g, 


dv/dt = -k/m v2 + g’ (1) 


dv/dt = -k/m(v2 — mg’/k) 
Letting dv/dt = 0 to determine the drone’s terminal speed, 
0 = -k/m(vz — mg’/k) 
vr = V(mg’/k) (2) 
Returning to (1), which is separable and integrable, 
dv/(v2 — mg’/k) = -k/m dt 
Using (2), 
k/m = g'/vt 


dv/(v2 — v4) = -k/m dt with v(0) = 0 


Vv t 
| aviv? = v3) = - e709 | at (3) 
0 0 


Using the method of partial fractions for the integrand on the left side of (3), 
1/(v2- va) = A/(v - vr) + B/(v + vr) (4) 
After a bit of algebra, 
Ov+1=(A+B)v+ (A-B)vr 
Comparing coefficients of like terms, 
(A+B)=0orB=-A 
and 
(A - B)vr = 1 = (A+A)vr=1 or A= 1/(2vr) 


Substituting the expressions for A and B into (4), 
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1/(v? —vp) = 1/(2vr)[1/(v - vr) - 1/(v + vr) (5) 


Substituting (5) into the integrand on the left side of (3), 


1/(2v7) i es - lv = 83 f a (6) 


(v— vr) 


Integrating (6) by inspection, and noting that v < vz, 
1/(2vr)[In((v - vr)/-vr) - In((v + vr)/vr)] = -g’t/vy 
Using properties of logarithms, 
In[-(v - vr)/(v + vr)] = -29’t/vr 
-(v- vr) /(V + vr) = e*(-2g't/vr) 
After some algebra, 
V =vr[1- e*(-2g’t/vr)]/[1 + e*(-2g’t/vr) | 
Multiplying the numerator and denominator by e’(g’t/vr)] (a cute ‘trick’), 
v= vile*(g't/vr) — e*(-g't/vr)]/[ e*(g't/vr) + e*(-g't/vr)] 
Noting that 
2sinh(g’t/vr) = e*(g’t/vr) — e*(-g’t/vr) 
and 
2cosh(g’t/vr) = e*(g’t/vr) + e4(-g’t/vr), 
then 


v = vrsinh(g’t/vr)/cosh(g’t/vr) 


v = vrtanh(g’t/vr) (7) 
Assuming g = 0.62 m/s?, 8 = 23°, m = 4.00 kg, k = 0.08 kg/m, and p = 0.10, 
g’ = [sin(23°) — 0.10cos(23°)|0.62 m/s? = 0.185 m/s? 
and 
vr= V(4 kg x 0.185 m/s2/0.08 kg/m) = 3.04 m/s 
Rewriting (1), using the above parameters, and omitting units for clarity, 
dv/dt = -0.020v2 + 0.185 (8) 
Well, you were warned that the math would be a bit lengthy! © 
3 Computer setup 
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Figure 1: Computer setup for navigating the slopes 


Modified Euler method 
(hand-held programmable calculator) 


PROGRAM:PLUTOV 
:ClrHome:ClrDraw 
’"DV/DT=-0.020V2 +0.185” 
“INITIAL CONDITION:” 
:0-T:0-V 

:’STEP SIZE:":1>H 


‘Fix 2 
:Lbl 1 


‘If T>25:Then 

:Goto 2:Else 

:Disp {T,V} 
:-0.020V2+0.185-F 
:T+HOT 
:V+FH>W 
:-0.020W2+0.185-S 
((F+S)/2A 
:V+AH-V 
:Pause:Gotol 

:Lb] 2:End 


4 Results 





Time 


Analog Computation Solution 
estimated from oscilloscope 


Analytical Solution 


Numerical Method Solution 
from hand-held calculator 
































t(s) | v(Gm/s) v(im/s) v(m/s) 
00.0 | 0.0 0.00 0.00 
10.0 | 1.6 1.65 1.65 
20.0 | 2.5 2.55 2.55 
30.0 | 2.8 2.89 2.89 
40.0 | 2.9 2.99 2.99 
50.0 | 3.0 3.03 3.03 
>50_ | terminal speed = 3.0 3.04 3.04 








Table 1: Solution Comparisons 
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Figure 2: Drone speed as a function of time* 
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Figure 3: Drone terminal speed (t > 50 s, not quite infinity!) via digital multimeter 


*For this application note, the display was produced during a single run by a 
differential equation analog computer prototype using discrete components 
with tolerances between 1% and 10%. 
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